1 



o 
o 

(N 



C/3 ' 

O ' 

• i-H ' 

-f— > ■ 

Or 

q ; 

o ■ 
^ ■ 



> ■ 

m : 
o 

m ■ 
o ■ 
m ' 
O . 

o : 
• i— i 

^ : 

Ph. 



I. INTRODUCTION 



Wave propagation in periodic layered media is very 
similar to the motion of electrons in crystallyne solids 
(hence, they are also known as photonic crystals). Thus, 
it is hardly surprising that some well-established con- 
cepts of solid-state physics (such as Bloch waves, Bril- 
louin zones, and forbidden bands) have been incorporated 
in the doctrine of photonic crystals 

The standard approach 2] identifies the periodic struc- 
ture with a one-dimensional lattice that is invariant un- 
der lattice translations, which requires the system to be 
strictly infinite. Of course, this is unattainable in prac- 
tice. In solid-state physics one goes around this difficulty 
by considering that the one-dimensional lattice closes 
into a ring; i.e., by imposing periodic boundary condi- 
tions. While this seems to be reasonable for a lattice of 
atoms, is obviously unfeasible for periodically stratified 
media. 

In spite of this fact, only under the assumption of an 
endlessly repeating stack of dielectric slabs this standard 
treatment gets an eigenvalue equation for the Bloch fac- 
tor and then one finds 0, EL El El 0] that strong reflection 
(stop bands) occurs when the trace of the 2x2 matrix 
describing the basic period exceeds 2 in magnitude. 

In practice not very many periods (say, ten) are needed 
to have a stop band. Therefore, one is led to consider 
stacks of N periods, which are often appropriately called 
"finite" periodic structures @. The matrix describing 
these structures can be computed as the Nth power of the 
basic period in terms of Chebyshev polynomials 0, . 
This gives an exact expression for the reflectance, but so 
complicated that it is difficult to get, at first glance, any 
feeling about the performance of the system. Moreover, 
the distinction between allowed and forbidden bands is 
obtained only numerically. 

The aim of this paper is to provide an alternative set- 
ting for dealing with these finite periodic structures that 
stresses the different scaling laws appearing in the re- 
flectance. To this end we resort to the Iwasawa decom- 
position to represent the action of any multilayer as the 
product of three fundamental matrices of simple interpre- 
tation. We show that these matrices present quite dis- 
tinct reflectances. The trace of the basic period allows us 
to classify periodic structures in three types with prop- 
erties closely related to one (and only one) of the three 
fundamental matrices. Finally, we study the reflectance 
associated to each one of these types. 

The advantage of this approach is, in our view, twofold: 
it naturally emphasizes the different laws for the for- 
bidden bands, allowed bands, and band edges and, at 
the same time, provides universal expressions remarkably 
simple for the corresponding reflectances. 





FIG. 1: Wave vectors of the input [Ea and E^ '] and the 



output [El > and Ey'] fields in a multilayer sandwiched be- 
tween two identical semi-infinite ambient and substrate me- 
dia. 



II. THE ORIGIN OF THE BANDGAP IN 
INFINITE PERIODIC STRUCTURES 



To keep the discussion as self-contained as possible, we 
first briefly summarize the essential ingredients of mul- 
tilayer optics that we shall need for our purposes. The 
basic period of the structure is illustrated in Fig. 1 and 
consists of a stack of 1, . . . , j, . . . , m plane-parallel layers 
sandwiched between two semi-infinite ambient (a) and 
substrate (s) media, which we shall assume to be iden- 
tical, since this is the common experimental case. Here- 
after all the media are supposed to be lossless, homoge- 
neous, and isotropic. 

We consider an incident monochromatic linearly polar- 
ized plane wave from the ambient that makes an angle 9q 
with the normal to the first interface and has amplitude 
E a . The electric field is either in the plane of inci- 
dence (p polarization) or perpendicular to the plane of 
incidence (s polarization). We consider as well another 
plane wave of the same frequency and polarization, and 



?(-) 



incident from the substrate at the 



with amplitude Ei 
same angle Oq. 

As a result of multiple reflections in all the inter- 
faces fffj ]. we have a backward-traveling plane wave in 

the ambient, denoted Ea , and a forward-traveling plane 



wave in the substrate, denoted E 



(+) 



If we take the field 



2 



amplitudes as a vector of the form 

25<+> 
C 1 £?(") 



(1) 



which applies to both ambient and substrate media, then 
the amplitudes at each side of the multilayer are related 



by the 2x2 complex matrix I 
matrix, in the form 



we shall call the transfer 



E a = M as E s . (2) 
The matrix M as can be shown to be of the form [l2l[T^ | 



Mas 



1/Ias Ras/Tas 
Ras/T as 1/T as 



a (3 
(3* a* 



where the complex numbers 



R as = \R as \ exp(ip), T as = \T as \ cxp(ir), 



(3) 



(4) 



are, respectively, the overall reflection and transmission 
coefficients for a wave incident from the ambient. Note 
that we have det M as = +1, which is equivalent to 
|-Ras| 2 + \T a s\ 2 — 1) an d then the set of lossless multilayer 
matrices reduces to the group SU(1,1), whose elements 
depend on three independent real parameters [Til Il5| . 

The identity matrix corresponds to T as — 1 and R as = 
0, so it represents an antireflection system without trans- 
mission phase shift. The matrix that describes the over- 
all system obtained by putting two multilayers together 
is the product of the matrices representing each one of 
them, taken in the appropriate order. So, two multilay- 
ers, which are inverse, when composed give an antireflec- 
tion system. 

When we repeat the basic period M as the fields at the 
nth cell are given by 



E„ 



(5) 



In an infinite structure the Floquet theorem ensures that 
these cells are equivalent and so the two vectors in Eq. (JHJ 
are proportional 0, H] 



E n 



XE„ 



(6) 



and therefore the Bloch factor A is determined as the 
eigenvalue of the basic-period transfer matrix M as : 



which gives 



(M as - X)E n = 0, 



A 2 -[Tr(M as )]A+l = 0. 



(7) 



(8) 



If [Tr(M QS )] 2 > 4, the solution grows or decays exponen- 
tially and no propagating wave is possible (stop band). 
On the contrary, Bloch waves appear in this algebraic 
approach when [Tr(M as )] 2 < 4 and then the matrix M as 
has two eigenvalues (allowed band). The band edges 



are obtained by setting [Tr(M as )] 2 = 4. Such a peri- 
odic structure exhibits resonance reflection very much 
like the diffraction of X rays by crystal lattice planes and, 
in consequence, it is usually called a Bragg reflector 0. 
By properly designing the basic period it is possible to 
achieve extremely high reflectance for some selected spec- 
tral region ^(| . 



III. FINITE PERIODIC STRUCTURES AND 
IWASAWA DECOMPOSITION 

As anticipated in the Introduction, an infinite multi- 
layer is unrealistic in practice, and one should consider 
iV-period finite structures, for which the previous stan- 
dard approach fails [13, Ell • 

The overall transfer matrix for a TV-period struc- 
ture is equal to (M as ) N . Explicit but complicated for- 
mulas for the corresponding matrix elements exist pj.llC|. 
but they do not give analytic expressions for the bandgap 
and are difficult to analyze. 

We have recently advocated us ing the Iwasawa decom- 
position to study layered media since it provides a 
remarkable factorization of the matrix representing any 
multilayer (no matter how complicated it could be) as 
the product of three matrices of simple interpretation. 
The decomposition can be stated as 



where 



m = 

HO = 
N 0) = 



\ as = K(0)A(ON(i/), 



exp(i^/2) 

exp(-i(j)/2) 

cosh(£/2) isinh(£/2)\ 
-isinh(£/2) cosh(£/2) ) ' 

l-iv/2 v/2 
v/2 l + iv/2 



(9) 



(10) 



The parameters <f>, £, and v are given in terms of the 
elements of the multilayer matrix by 



4/2 = arg(a + %0) , 
i/2 - ln(l/|a + i/3|) , 
v/2 = Re(aP*)/\a + i(3\ 2 



(11) 



and their ranges are and — 2ir < <j) < 2ir. Now, 

we can interpret the physical action of the matrices ap- 
pearing in Eq. Q. K(<f>) represents the free propagation 
of the fields E in the ambient medium through an optical 
phase thickness of <f>/2, which reduces to a mere shift of 
the origin of phases. Alternatively, this can be seen as an 
antireflection system. The second matrix A(£) represents 
a symmetric system with real transmission coefficient 
Ta = sech(£/2) and reflection phase shift /5a = ±7r/2. 
Finally, the third matrix, N(^), represents a system hav- 
ing Tim = cos(tn) exp(ir N ) and i?N = sin(r N ) exp(ir N ), 



3 



with tan(ri\i) = v/2. In Ref. [l^ we have provided sim- 
ple examples of realistic systems that accomplish these 
requirements. 

For our purposes here, the essential point is that the 
reflectance 9\ — \R\ 2 associated to each one of these ma- 
trices is 



= 0, 

m A = tanh 2 (£/2), 

5K N = (z,/2) 2 /[l + (^/2) 2 ] 



(12) 



While 91k is identically zero, 9^a and 9^n tend to unity 
when £ and v, respectively, increase. However, they have 
distinct growth: 91a goes to unity exponentially, while 
£Hm goes as 0(v~ 2 ). 

The Iwasawa decomposition proves also to be a pow- 
erful tool for the classification of layered media. Indeed, 
after Ref. j2£j one is led to introduce the following cri- 
terion: a matrix is of type K when [Tr(M as )] 2 < 4, is 
of type A when [Tr(M QS )] 2 > 4, and finally is of type N 
when [Tr(M as )] 2 = 4. Although this trace criterion has 
an elegant geometrical interpretation [l3T | and coincides 
with the one giving the stop bands in Eq. JSJl, let us re- 
mark that if a multilayer has a transfer matrix M as of 
type K, A, or one can always find a family of matrices 
C [also in SU(1,1)] such that 



= C 



(13) 



is just a matrix K ((f)), A(£), or N(f), respectively. 
This conjugation by C preserves the trace: Tr(M as ) = 
Tr(M as ). Moreover, Eq. (|13f) allows us to recast the mul- 
tilayer action in Eq. J5J) as 



E a — M as E s 
where the new field vectors are 
E = CE. 



(14) 



(15) 



In other words, the matrix C gives a new vector basis 
such that the basic period of the system, when viewed in 
such a basis, presents a reflectance exactly of the form in 
Eq. (|12|l . This provides a full characterization of the re- 
flectance of any periodic system, as we shall demonstrate 
in next Section. 



IV. CHARACTERIZING THE THREE 
FUNDAMENTAL BEHAVIORS OF THE 
REFLECTANCE 



where \i represents the appropriate parameter cf>, £, or v. 
For a iV-period system the overall transfer matrix Mis 
is 

MW - C- 1 [M as (fi)] N C = (T 1 M as (N(x) C, (17) 

which does not depend on the explicit form of the basic 
period. 

After Eq. (|17|) , one must expect three universal behav- 
iors of the reflectance according the basic-period transfer 
matrix is of the type K, A or N. We shall work in what 
follows the detailed structure of these three basic laws. 

Because the stop bands are given by the condition 
[Tr(M as )] 2 > 4, we first consider the case when M as is 
of type A. Then, Eq. JTSJ gives 



C- 1 A(x) C, 



where we have denoted 



Re(a) = -Tr (I 



= cosh(x) > 1, 



(18) 



(19) 



because we are taking into account only positive values 
of Re(a), since negative values can be treated much in 
same way. If we put 



C 



Cl c 2 

one solution of Eq. (|18|1 is 

Ci = F(J3* + i sinhx), c 2 
where the factor F is 

F = 



(20) 



-iF1m(a), (21) 



y/2 sinh x[sinh \ I m (/5)] 



(22) 



Carrying out the matrix multiplications in Eq. i|17|) it 
is easy to obtain the reflectance of the A^-period system 
as 



(AO 



[sinh(x)/ smh(Nx)} 2 



(23) 



This is an exact expression for any value of JV. As JV 
grows, %Vfc approaches unity exponentially, as one could 
expect from a band stop. 

The band edges are determined by the limit condition 
[Tr(M as )] 2 = 4; that is, when M os is of type N. A calcu- 
lation very similar to the previous one shows that now 



Since the matrix M as associated to the basic period 
belongs to one of the subgroups K ((f)), A(£), or U(u) of 
the Iwasawa decomposition, and all these subgroups are, 
in our special case, Abelian and uniparametric, we have 
that 



M QS (Mi)M a;s (^ 2 ) = M aa (/xi + fx 2 ) 



(16) 



m 



(N) _ 
N ~ 



(1/A0 2 ' 



(24) 



with a typical behavior ~ 1 — 0(N~ Z ) that is uni- 

versal in the physics of reflection. The general results 
(|23|l and l|24l) have been obtained in a different frame- 
work by Yeh and Lekner Q . 
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Finally, in the allowed bands we have [Tr(M as )] 2 < 4; 
M as is of type K, and 



{N ) _ Q 2 - 2Qcos(2jVO) 
~ 1 + Q 2 - 2Qcos(2AT7) : 



where 



and the phase O is determined by 



(Re a) + iy/l - (Re a) 



(25) 



(26) 



(27) 



Now the reflectance oscillates with N between the values 
(Q 2 - 2Q)/(Q - l) 2 and (Q 2 + 2Q)/(Q + l) 2 . 

Equations (|23ll , (|24H and 125() are the three fundamen- 
tal reflectances we were looking for. In our opinion, they 
have the virtue of simplicity. In practice, the basic pe- 
riod usually consists of two homogeneous layers and the 
expressions become simpler 0. We emphasize that the 



scaling laws expressed by these equations are universal 
features of the reflection in finite periodic systems. 



V. CONCLUSIONS 

We have used a trace criterion to provide in a direct 
way three fundamental forms for the reflectance that are 
universal for any multilayer. When the system is finite 
periodic this result has allowed us to obtain in a direct 
way explicit expressions for the reflectance of the overall 
system. 

This approach has physical meaning on its own, irre- 
spective whether the system is periodic or not. Moreover, 
it leads naturally to a further analysis of finite periodic 
bandgap structures that does not need the strict require- 
ment of infinite periodicity, as it is assumed in the stan- 
dard theory. 

We acknowledge illuminating discussions with Jose F. 
Carinena. 
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